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Abstract 


Given the complicated mechanisms associated with the operation of the solid oxide fuel cell (SOFC), mathematical models developed to capture 
the SOFC’s dynamic characteristics often result in high order and complex dynamics that make the model unsuitable for control design and analysis. 
In this paper, the minimum Gibbs free energy (MGFE) method is exploited to simplify the calculation of the mass balance dynamics of the fuel 
flow in the SOFC in an effort to develop a control oriented model that achieves appropriate trade-off between model accuracy and simplicity. 
The simplified model is compared with a baseline model where dynamic governing equations are derived for the mass balance of each of the six 
gas species in the fuel flow. The implications of the MGFE method on the modeling performance are investigated through numerical simulations 
and frequency domain analysis. The MGFE method leads to the elimination of 5 out of 11 states in the baseline model, thereby resulting in a 
significantly lower order model. Critical parameters that may influence the accuracy of the simplified model are also identified. The study concludes 
that the accuracy of the reduced-order model is acceptable at normal conditions for both steady-state and transient operations. Noticeable model 
errors are only observed in steady-state responses at extreme conditions, where the fuel inlet has a high fraction of CH, and the fuel cell has low 
fuel utilization, and in transient when the inputs change at high frequency. 


© 2006 Elsevier B.V. All rights reserved. 
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1. Introduction 


Given its high efficiency, low emissions and flexible fueling 
strategies, the solid oxide fuel cell (SOFC) offers a promising 
alternative energy solution for a wide range of applications. As 
the SOFC technology keeps evolving rapidly, dynamic mathe- 
matical models become critical tools that facilitate the design, 
analysis and optimization of SOFC systems. Intensified research 
efforts have been made on the SOFC model development in 
recent years, and several dynamic models have been reported 
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[1-8]. Through model-based simulations, potential operation 
issues, such as the slow load following and large overshoots 
of the temperature and temperature gradient, have been iden- 
tified in the transient response of the SOFC [8-11]. Feedback 
control strategies are being explored to improve the SOFC’s 
performance [10,11]. 

Compared to the substantial progress achieved in the mod- 
eling, material and manufacturing of the SOFC, however, there 
has been much less results reported on the control design for 
SOFC systems in the literature. One important reason is prob- 
ably because the dynamic models are often too complex for 
applying the model-based control methodologies. In order to 
capture the inherently complicated phenomena associated with 
the SOFC, the governing equations based on the electrochemi- 
cal, thermodynamic and gas flow principles will easily drive the 
model to prohibitively high complexity and order. For instance, 
a feedback controller is designed in [11] to speed up the load 
following response of an SOFC system, which has 196 states in 
the original model. Model linearization and balanced truncation 
are adopted to reduce the original model down to a seven-state 
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Nomenclature 

A area of the electrochemical reaction (m?) 

Ac cross area of the gas channel (m?) 

Cp heat capacity of solid layers (J K7! kg~!) 

Cys heat capacity of gaseous species s (J K7! mol7!) 

Cs molar concentration of species s (mol m~?) 

d channel height (m) 

es specific internal energy of species s (J mol~!) 

F Faraday constant (C mol`!) 

hs(T) specific enthalpy of species s at temperature T 
(Jmol!) 

i current density (A m~?) 


kf sol» Kao! heat transfer coefficients between bulk flows 
and their solid walls (J K7! m~? s7!) 

l cell length (m) 

M; molar mass of species s (kg mol” !) 

Nin,s, Nout, inlet and outlet molar rate of species s, respec- 
tively (mol s7!) 


Ps partial pressure of gas species s (Pa) 
P total pressure (Pa) 
Po downstream pressure of SOFC (Pa) 


qin, Gout inlet and outlet enthalpy flux of gas flows, respec- 
tively (J s7! m7?) 


Tk rate of reaction k (mol s~! m7?) 

R universal ideal gas constant (J mol`! K7!) 
T temperature (K) 

Uout outlet flow velocity (m s7!) 


U, Uocv operating voltage and open circuit voltage of the 
fuel cell (V) 


V volume of gas channel (m?) 

Wout outlet mass flow rate (kg s7!) 

Xs molar fraction of species s 

Greek symbols 

a orifice constant of gas channels (ms) 

Ver fuel utilization ratio 

Yor pre-reforming ratio 

YO./C Oxygen-to-carbon ratio 

n potential loss (V) 

s/c molar ratio of s/C of the bulk flow in the fuel 
channel 

Vs,k stoichiometric coefficient of species s in reaction 
k 

p density (kg m7?) 

T solid layer thickness (m) 

Subscripts 

a air flow 

act activation 

con concentration 

f fuel flow 

I interconnector 

ohm ohmic 

Ox oxidation reaction 


PEN PEN structure 


Red reduction reaction 

Sa species in the air flow 

Sf species in the fuel flow 
sol solid structure in SOFC 
SR steam reforming reaction 


WGS water gas shift reaction 


linear one, before the linear quadratic Gaussian (LQG) method 
can be applied to derive feedback control laws. The disadvan- 
tages of this approach, however, include the limited applicable 
operating range of the resulting controller as well as the difficulty 
of relating the states in the reduced-order model to physical vari- 
ables. In order to make the model amenable for existing control 
design platforms, control-oriented models, which preserve the 
key dynamic characteristics of the system over a wide operating 
range while remain low order, are in high demand to make the 
model-based control development possible. 

One big challenge in simplifying the SOFC model attributes 
to the complicated composition and reactions in the bulk flow 
along the fuel channel. For example, while some models [2,4,6] 
assume that the fuel inlet to the SOFC contains only H? and H20, 
others [12,13] allow the carbon monoxide in the reformate enter- 
ing the SOFC. In fact, CO is considered as fuel for the SOFC, 
either directly participating in the electrochemical reaction or 
producing hydrogen through the water gas shift (WGS) reac- 
tion. Direct internal reforming (DIR) is also proposed to allow 
hydrocarbons to be reformed directly inside the fuel channel 
of the SOFC [12]. Whilst these fuel flexibilities bring bene- 
fits of reducing the size of the external reformer and increasing 
the efficiency of the system, the complex gas compositions and 
chemical reactions involved lead to complicated calculations of 
the mass and energy balances in the SOFC. In [8], adynamic gov- 
erning equation is derived for the mass balance of each of the five 
species in the fuel flow, accounting for 5 out of 12 states in each 
discretization unit used in the finite volume method. In [1,5], 
the mass and energy balance in the gas flows are approximated 
in a quasi-static manner, resulting in a set of coupled non-linear 
algebraic equations that usually require extensive computational 
resources to solve. Although the SOFC is often modeled as a 
distributed parameter system, it should be noted that the issue 
mentioned above still exists even when the lumped parameter 
assumption is adopted. 

The minimum Gibbs free energy (MGFE) method is an 
approach to obtain the equilibrium composition of chemical 
reactions, by minimizing the Gibbs free energy of the mix- 
ture [14]. In this paper, we use the MGFE method to simplify 
the calculation of the gas composition in the fuel channel and 
investigate the implication of this model reduction approach. 
For simplicity, the SOFC is considered as a lumped parameter 
system. A baseline model is first developed, where dynamic gov- 
erning equations of mass balance are derived for each species in 
the fuel flow. A reduced-order model is then developed using the 
MGFE approach and compared to the baseline model for both 
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Interconnector 
CH4, CO», CO,| CHa +H,O — 3H, +CO (steam reforming) CHa, CO2, CO 
H20, H2, N2 |CO +H,O—> H, + CO (Water gas shift) H20, Ho, No 
Fuel flow H2 HO Fuel channel ; 
e 
O2, N2 
> Oz Air channel c> 
Air flow one 
Interconnector 


Fig. 1. Operating principle of co-flow planar SOFCs. 


steady state and dynamic performance. Critical parameters that 
may affect the validity of the MGFE model are also identified. 
The rest of the paper is organized as follows: the baseline 
model is described in Section 2 and a reduced-order model using 
the MGFE method is derived in Section 3. Based on simula- 
tion results of both steady-state and dynamic responses, these 
two models are then compared in Section 4, in order to investi- 
gate in depth the benefits and potential pitfalls of the simplified 
model. Two fuel processing methods and different operating 
conditions are considered in analyzing the impacts of using the 
MGFE method on the accuracy of the simplified model under 
different scenarios. Both steady-state and transient operations 
are considered in the analysis. The conclusions are given finally. 


2. Modeling of planar SOFCs 


The operating principle of the planar SOFC considered in this 
paper is illustrated in Fig. 1. In order to make the model com- 
patible with different system configurations, such as different 
types of fuel processors and anode recirculation schemes, the 
fuel inlet to the fuel cell can be any combination of CH4, CO2, 
CO, H20, Hə and Na, and direct internal reforming is included 
in the model. Table 1 lists all the reactions to be considered in 
the model. 

To focus on the implications of the MGFE method on model 
performance, the SOFC under consideration is treated as a 
lumped parameter system by neglecting the spatial distributions 
of variables, such as the current density, gas compositions, tem- 
peratures and pressures, along the direction of the bulk flows in 
this paper [6]. These variables are assumed to be homogeneous 
in the fuel cell, and the dynamic governing equations are derived 
by applying the electrochemical, thermal dynamic and gas flow 
principles. It should be noted that the proposed approach is also 
applicable to situations where the spatial distributions of these 
variables are considered. In the latter case, the system is often 
divided into smaller units and the governing equations can be 
applied to each discretization unit of the whole cell. 


Table 1 

Reactions considered in the model 

Location Reaction Expression 

Fuel channel SR CH4 + H20 — CO + 3H2 
WGS CO + H20 = CO2 + H2 

Anode Ox Hə + O% > H20 + 2e7 

Cathode Red 0.502 + 2e7 > 0%- 


For easier cross reference, assumptions used to develop the 
baseline model are summarized as follows: 


(1) Current is produced only by oxidation of H2, and CO only 
reacts through the WGS reaction. 

(2) Ideal gas flows for both fuel and air channels. 

(3) Constant Nusselt number. 

(4) Adiabatic boundaries for the cell. 

(5) Outlet gas flows have the same composition and temperature 
as inside the gas channels. 


2.1. Electrochemical sub-model 


The operating voltage of the cell can be calculated as follows: 


U = Uocv — (ohm + Nact + Neon); (1) 


where the last three terms account for various potential losses. 
Uocv is the open circuit voltage determined by the Nernst equa- 
tion: 
RTpEn DHO 
Uocv = Eo - —;p In f aE: 
Pp H2 P Oo 


(2) 


with Eo = 1.2723 — 2.7645 x 10~4Tpen [15], where Tpgn is 
the temperature of the PEN structure, and pupo, pu, and po, 
are the partial pressures of H20, Hz and Oz, respectively. 

The activation loss, nact, is due to the energy barriers to be 
overcome in order for the electrochemical reaction to occur, 
and can be characterized by the Butler-Volmer equation with 
the transfer coefficient of 0.5 [16]. The concentration loss, Neon, 
reflects the overpotential due to the species diffusion between 
the reaction site and the bulk flow in gas channels, and can be 
calculated by the approach used in [7]. Nohm is the ohmic loss 
due to the electrical and ionic resistance along the path of the 
current in the fuel cell. 


2.2. Mass balance sub-model 


In this sub-model, dynamic equations are derived to calcu- 
late species concentrations of the bulk flows in the fuel and air 
channels. The mass balance dynamics in fuel channel can be 
described as follows: 


; 1 1 
Cs = (Nin, s¢ = Nours) qr + 5 Vsg,k Tk 
f  ke{SR,WGS,Ox} f 
sp € {(CH4, CO2, CO, H20, H2, No}, (3) 


where C;, is the molar concentration of species s¢ in the fuel 
flow, Nin,s, and Nout,s, the inlet and outlet molar flow rates of 
species ss, respectively, Vp the volume of the fuel channel, vsp, 
the stoichiometric coefficient of species sp in reaction k, rg the 
rate of reaction k, and df the height of the fuel channel. 

According to Faraday’s law, the rates of Ox and Red reactions 
are related to the current density as follows: 

i 


TOx = Red = 55 (4) 
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Although different expressions are proposed to calculate the 
reaction rate of the SR reaction, the following formula is often 
used to represent the typical kinetics of the internal SR reaction 
in the SOFC [1,7,15-17]: 


FSR = ksR PCH, XP (-=) , (5) 
RT; 
with ksr = 0.04274 mol (s m? Pa)~! 
mol™!. 
The WGS reaction is usually considered a very fast one and 
is assumed to be at its equilibrium [1]. The formula given in [7] 
is used here to account for this effect: 


PCO PH) ) 
PCO PH:0 Keg,wes / ` 


and Esr = 82, 000J 


rwos = kwas Pco (1 (6) 
where kwos = 0.01 in this model and Keg,wes the equilibrium 
constant with Keg, wos = exp(4276/T; — 3.961) [17]. 

While Nin,sẹ is specified by the fuel inlet condition, Nout, s 18 
determined by the following relation of the fuel flow continuity: 


Nout, s¢ = Uout,fCstAc,f» (7) 


where uout, is the speed of the outlet fuel flow and Ac + the cross 
area of the fuel channel. Given the small pressure drop across 
the fuel cell, the linear orifice relation is adopted to calculate 
Uoutf as follows: 


Wout, 
fee (8) 
E Ac,f Žas Ms;Cs¢ 
Wout,t = &fCPr — Po), (9) 


where Wout is the outlet mass flow rate of the fuel flow, Ms; the 
molar mass of s¢, œ¢ the orifice constant of the fuel channel, Pr 
the pressure of the bulk flow in the fuel channel and Pp is the 
downstream pressure of the SOFC. Plugging Eqs. (8) and (9) 
into Eq. (7), we have: 
ap(Pe — Po) 
a Ms; Cse 
The partial pressures of species and total pressure of the fuel 
flow can be determined as follows: 


Psg = RT;Cs,, (11) 


P=) Poss (12) 


Se 


Cy. (10) 


Nout,s¢ = 


where 7; is the temperature of the fuel flow, which is solved in 
the energy balance dynamics. 
Similarly, the mass balance in the air flow follows: 


: 1 1 
Cs, = (Nin, s, = Nout,sa) > + Vs,,Red/Red—, Sa € {O2, No}, 
Va da 


(13) 
with 


i RRO, (14) 


a Ps (15) 
Sa 


&a(Pa — Po) 
= —— 16 
Uout,a Aca = M,,Cs, (16) 
aa(Pa — Po) 
Nout,sq = Uout,a Cs, Aca = tc (17) 


E Ze Ms, Csa ™ 
2.3. Energy balance sub-model 


As in [1,5], the cell is divided into three temperature layers, 
i.e., the fuel/air bulk flows (7¢/7T,) and solid structure (Tsoi), 
where the solid structure includes both the PEN and intercon- 
nectors, i.e., TPen = Ti = Tsoi. The temperatures of these layers 
are calculated by solving the dynamic equations of the energy 
balance in each layer. 

The energy balance dynamics in the fuel flow can be 
expressed as follows: 


d 1 1 
dt 2 Cstêsp = (qin,f E dout t) 7 F 2kfsol(Tsol _ a 


1 
+ roxlhn,o(Tpen) — hen TT (18) 


where es; is the specific internal energy of species sp, the first 
term on the right hand side of Eq. (18) is due to the enthalpy flux 
of the bulk flow, the second term accounts for the convective 
heat exchange between the fuel flow and its surrounding solid 
layers. The heat transfer coefficient can be obtained by assuming 
constant Nusselt number of 4 [16]. The last term of Eq. (18) is 
caused by the enthalpy flux due to the Ox reaction at the anode. 
By the relation for the ideal gas flow, és, = Asp — (Ds-/Cs¢), 
where Asẹ is the specific enthalpy of species ss, we can obtain: 


a Sos (TRC: + (qi i 
f= a” oe - spl Lf. f) sp din,f — Yout,f 7 


1 1 
+ 2k¢ so1(Tso1— Th) g trollo) — hy, (74)] P | ; 


(19) 


Similarly, for the air flow, we have: 


i= DETA {- 2 Mie) = RT IG + (Gina—out,a) 


1 1 1 
aa 2ka,sol(Tsot — T,)— — 0.5rRea ho, (Ta)— : (20) 
l da da 


The inlet enthalpy flux, ging and gina, are dependent upon 
the fuel and air inlet conditions, respectively, while the outlet 
enthalpy flux of the bulk flows can be calculated by: 


dout,f = Uout,f 5 Csehse(Te), (21) 


SE 


Gouta = touta X Csahs, (Ta). (22) 


Sa 


H. Xi et al. / Journal of Power Sources 165 (2007) 253-266 257 


The energy balance dynamics in the solid structure can be 
derived as follows: 


; 1 
Tsol = 


PPENC p,PENTPEN + PIC p,ITI 
x { —2k¢.so1(Tsol — T- 2ka,sol(Tso1 — Ta) 
+ rox [hi (T})+0.5h0,(Ta)—hH,0(Tso1)] — iU}. 23) 


3. Model reduction for planar SOFCs using MGFE 
method 


Itis noted that the dynamic SOFC model developed in Section 
2 has 11 states. Among them, six states, which represent the 
concentrations of the six species in the fuel flow, come from the 
mass balance dynamics of the bulk flow in the fuel channel. In 
this section, a reduced-order model is derived using the MGFE 
method to simplify the calculation of the composition and mass 
balance in the fuel flow, leading to a control-oriented SOFC 
model with a substantially lower order. 


3.1. Composition calculation based on MGFE method 


By the MGFE approach, the equilibrium composition of the 
chemical reactions can be obtained by minimizing the Gibbs free 
energy of the mixture of reactants and products [14]. Assume 
we consider a generic reaction as follows: 


aA+bB=dD+eE, (24) 


and have nl, s € {A, B, D, E} moles of species s at the begin- 
ning. Denote n¥ as the amount of species s at the equilibrium 
of the above reaction at temperature T. Following the MGFE 
method, we have: 


Sie = min $ nsG', s €{A, B, D, E} (25) 
S AY 
0 0 0 0 
subject to TATTA S B R TDT TD ETE (26) 
—a —b d e 


where GT is the Gibbs free energy of species s at temperature 
T. The constraint given in Eq. (26) follows the law of mass 
conservation. 

The MGFE method is particularly useful when multiple 
chemical reactions occur simultaneously, such as the SR and 
WGS reactions in the fuel channel in our case. Five species, 
namely CH4, CO2, CO, H20 and Hyp, are involved in these two 
reactions. Since the amount of the element of carbon (C) is 
invariant before and after the reactions considered in the fuel 
channel, the MGFE method for the SOFCs fuel flow can be 
formulated in terms of ratios of chemical species. Define: 


N. 
Bsa/c = i so € {C, H, O, CH4, CO2, CO, H20, H2}, (27) 
Cc 
where Ns is the amount (in moles) of so in the fuel channel. The 
normalization process of Eq. (27) is beneficial, as it will reduce 
the number of variables that have to be considered to determine 
the mixture’s composition. 


Denote cas jor 1 € {CH4, CO2, CO, H20, H2} as the value of 
Os, /c at the equilibrium of the SR and WGS reactions. By the 
MGFE method, given @/c, 90/c and the fuel flow temperature 
(Tr), Gy, jc satisfies: 


S onca = min S bacas, 
S1 s/c s1 
sı € {CH4, CO2, CO, H20, H3} (28) 


subject to @cHy/c + 9co,/c + 8co;c = 1, 
46cHu/C + 26H,0/c + 20H /c = On/c, 
26c0,/c + 8co/c + H,0/c = 9o/c; (29) 


where the constraints described in Eq. (29) reflect the 
mass conservation in terms of individual elements in the 
reactions. 

The fuel temperature, 7;, can be obtained from the energy 
balance dynamics of the fuel flow. Considering the small vol- 
ume of the fuel channel and neglecting the fuel accumulated 
therein, O4/c and 6o/c in the fuel flow can be calculated by 
the fuel inlet condition and current density. In fact, O4/c solely 
depends on the composition of the fuel inlet, while oxygen ions 
migrate from cathode to anode and enter the fuel flow through the 
electrochemical reactions. Define Ninc = Nin,cH, + Nin,co, + 
Nin,co. we have: 


_ 4Nin,cHy + 2Nin,H20 + 2Nin, H 


Onc i (30) 
Nin,C 
2Nin,co, + Nin,co + Nin,mo + U/2F) 
o/c = Rae . (31) 
1n, 


Given the fuel inlet composition and the current drawn from 
the SOFC, @yc and ĝo,;c can be calculated by Eqs. (30) and 
(31). This allows us to represent the MGFE results for A jc as 
functions of only three independent variables: 7;, Oy/c and 60/c. 

In this paper, the minimizers of Eq. (28), (A c are calculated 
off-line using a software named GASEQ [18], and the results 
are saved as look-up tables that can be used on-line in simula- 
tions. Examples of these look-up tables obtained are illustrated in 
Fig. 2. 

Once a jc are determined, the molar fraction of individual 
species in the fuel channel bulk flow, denoted as X;,, can be 
computed as follows: 


Nin,c%, jc 
Aa Fy 6 a+ NNy.in’ 
in,C ae s1/C F N2,in 
sı € (CH4, CO2, CO, H20, H2}, (32) 
N P 
Xn, 2 N2,in (33) 


Nin,C ee 0 /C T NN;,in 
3.2. Reduced-order SOFC model using MGFE method 


Instead of establishing a dynamic mass balance equation for 
each species in the fuel flow as been done in Section 2, we can 
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8 wo = 4,7 = 900 K 


0 He =8,T=900K 


o/c 


Owe =4,T=1100K 


—#— cue 


=m = O oHe 


8o/c 


Fig. 2. Examples of look-up tables obtained by the MGFE method. 


derive the following mass balance equation for the total mass of 
the bulk flow in the fuel channel: 
I 


me = Wins — Woutf + ap Mo» (34) 


where mg is the total mass accumulated in the fuel channel and 
the last term on the right hand side accounts for the migration 
of the oxygen ions from the cathode to anode. From Eq. (9), we 
have Wout, = &t( Pe — Po) where the total pressure of the bulk 
flow in the fuel channel can be calculated as follows in this case: 
me Š 

Pr Vi D X, M, RT;, (35) 
and consequently, the species partial pressures are Psp = PfX sp, 
where Xs; is determined by Eqs. (32) and (33) based on the 
MGFE method. 

In this paper, we focus on simplifying the mass balance in 
the fuel flow using the MGFE method as described above, and 
keep all other parts in the model given in Section 2 unchanged 
in order to investigate the implication of the MGFE approach. 
Further model simplification is possible. For example, equations 
similar to Eqs. (34) and (35) can also be derived for the air 
flow to remove one more state without having to use the MGFE 
approach. Furthermore, due to the small volumes of the gas 
channels, the mass balance and energy balance in the gas flows 
may be approximated in a quasi-static way to remove all the 
states in the mass and energy balance dynamics of the gas flows 
[1,5]. Given the scope of this paper, these model simplification 
approaches have not been applied. 

For the convenience of cross reference, the model developed 
in Section 2 is called the baseline model and the one in this 


section the MGFE model in the sequel. The MGFE method leads 
to the elimination of five states from the fuel flow mass balance 
dynamics of the baseline model. 


4. Model comparison 


The main difference between the MGFE model and the base- 
line model is in the composition calculation for the fuel bulk 
flow. For the MGFE model, the composition is assumed to be 
always at the equilibrium of the SR and WGS reactions and 
is determined using the MGFE method. The reaction kinetics 
are ignored in the MGFE model but are considered in the base- 
line model. The WGS reaction in the SOFC is considered fast 
[1,5,15,16], which has been taken into account in Eq. (6). The 
SR reaction as described in Eq. (5), however, is a relatively slow 
one and may take longer time to reach its equilibrium. There- 
fore, prediction errors may be expected for the reduced order 
model derived using the MGFE approach. Impacts of the MGFE 
assumption on modeling results are investigated in this section 
by comparing the MGFE model with the baseline model through 
simulations and analysis. 


4.1. Analysis of steady state performance 


Steady-state results obtained using the MGFE and base- 
line models are first compared, followed by the discussion on 
dynamic responses in the next subsection. 

Different types of fuel processors and operating conditions 
lead to different composition, particularly the fraction of CH4, 
in the syngas fed into the fuel channel of the SOFC, which could 


H. Xi et al. / Journal of Power Sources 165 (2007) 253-266 259 


Table 2 
Operating parameters of SOFC system 


Case number 


1 2 
Fuel processing method SR CPOX 
Current density (i) (A cm7?) 0.3, 0.6 
Fuel utilization ratio (yfr) 0.2, 0.4, 0.6, 0.8 
Fuel inlet temperature (K) 973 
Steam-to-carbon ratio 2 - 
Pre-reforming ratio (Ypr) 0.2, 0.5, 0.8 — 
Oxygen-to-carbon ratio (yo, /c) — 0.5, 0.6, 0.7 
Air excess ratio vi 
Air inlet temperature (K) 973 


dictate to what extent the underlying MGFE assumption is valid. 
Two cases as listed in Table 2 are considered in the following, 
where the fuel inlet to the SOFC are the resulting reformate using 
two different processing methods, one with steam reforming 
(SR) and the other using catalytical partial oxidation (CPOX). 
CH; is used as the fuel feedstock in both cases. 


4.1.1. SOFC with fuel processed by SR 
In the first case, we assume that CH4 is reformed in an external 
SR processor prior to the SOFC. The operating parameters of 


Current density=0.3A/cm? 


the system to be considered are listed under Case 1 in Table 2. 
For simplicity, the inlet temperatures of the fuel and air supplies 
to the SOFC, the steam-to-carbon ratio of the SR fuel processor 
and the air excess ratio are fixed. Different current densities 
(i), fuel utilization ratios (yfr) and pre-reforming ratio (Vpr) are 
considered in simulations, where yp; is defined as follows: 


amount of CH4 reformed in the processor 


Yp = (36) 


amount of CH4 supplied to the processor ` 
More CH4 will remain in the reformate at lower ypr. As the 
DIR is possible in the SOFC, partial reformation of CH4 in the 
reformer is considered and the remaining CH4 in the reformate 
reacts through the DIR in the fuel channel of the SOFC. 

Simulations are performed to obtain the steady-state 
responses of the SOFC at different operating conditions listed in 
Table 2. Fig. 3 compares the steady-state cell voltage and solid 
structure temperature (T7501) predicted by the baseline model 
and the MGFE model. Analyzing the simulation data, one can 
conclude that: 


(1) While the MGFE model produces significant errors at low 
fuel utilization ratio and low pre-reforming ratio, the two 
models match well at high yf or Ypr. Low fuel utiliza- 
tion ratios such as yfr = 20% lead to low SOFC efficiency 
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Fig. 3. Steady-state voltage and temperature of SOFC in Case 1. 
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Table 3 
Gas composition in fuel channel of SOFC in Case 1 (Xn, = 0) 
i (A cm~?) Ypr Ver Model Xcu, Xco, Xco Xpo Xp, 
0.8 Baseline 0.0013 0.1643 0.0348 0.6793 0.1202 
0.8 MGFE 0 0.1657 0.0343 0.6743 0.1257 
0.2 Baseline 0.0089 0.0787 0.1160 0.2967 0.4997 
0.3 MGFE 0.0005 0.0808 0.1189 0.2799 0.5199 
0.8 Baseline 0.0082 0.1697 0.0253 0.6923 0.1044 
0.2 MGFE 0 0.1700 0.0299 0.6700 0.1300 
0.2 Baseline 0.0636 0.0965 0.0654 0.3729 0.4017 
MGFE 0.0095 0.1021 0.0922 0.2744 0.5218 
0.8 Baseline 0.0022 0.1630 0.0357 0.6830 0.1161 
0.8 MGFE 0 0.1642 0.0358 0.6758 0.1242 
0.2 Baseline 0.0126 0.0757 0.1167 0.3060 0.4889 
0.6 MGFE 0.0003 0.0778 0.1220 0.2826 0.5172 
0.8 Baseline 0.0126 0.1688 0.0237 0.7049 0.0901 
0.2 MGFE 0 0.1683 0.0317 0.6717 0.1283 
0.2 Baseline 0.0785 0.0917 0.0612 0.4033 0.3653 
MGFE 0.0056 0.0976 0.0990 0.2721 0.5257 


and therefore are usually not expected in practical appli- 
cations. 

(2) Ts, obtained by the MGFE model is lower than that by the 
baseline model at all the operating conditions considered in 
Fig. 3, although the difference is insignificant at high yfr and 
Ypr- The model error exhibited in Tyo decreases as yfr OF Ypr 
increases, while it increases as i increases. 

(3) Except for yg, as high as 80%, the cell voltage predicted 
by the MGFE model is lower than that by the baseline 
model. For yfr < 0.6, the difference between the two models 
decreases as yfr increases. A decreased ypr or an increased 
i will lead to larger error in the steady-state cell voltage 
obtained by the MGFE model. 


In order to further identify the root cause of the modeling error 
introduced by the MGFE method, the molar fractions of the gas 
species in the fuel channel obtained by the baseline model and 
the MGFE model are compared in Table 3. One can see that, 
under the conditions considered, the Xcy, in the MGFE model 
is always lower than that in the baseline model, suggesting that 
the SR reaction be over-counted by the MGFE method, a direct 
result of the equilibrium assumption imposed by this approach. 

Due to the strongly endothermal nature of the SR reaction, 
the over-counted SR reaction in the MGFE model leads to 
the decreased temperature in the fuel cell as shown in Fig. 3. 
When Ypr decreases, the inlet fuel flow entering the SOFC has 
an increased fraction of CH4. Therefore, the different reaction 
kinetics adopted in different models have more significant influ- 
ence on the calculation of the fuel flow composition, leading to 
larger difference in XcyH, and consequently in Tsoi. This effect 
is exaggerated at lower yfr and higher i since the faster fuel 
flow involved reduces the time for CH4 to react through the SR 
reaction in the fuel channel. 

The cell voltage is affected by the MGFE approach primarily 
in two competing ways: on one hand, the decreased temperature 
in the SOFC will lead to reduced voltage; on the other hand, how- 
ever, the high partial pressure of H» resulted from the increased 


molar fraction of H2 in the anode due to the over-counting of 
the SR reaction will increase the cell voltage. The second fac- 
tor is dominated by the first one in most operating conditions 
considered in Fig. 3. At ye = 0.8, the second factor, however, 
becomes dominant where the low pressure of the hydrogen leads 
to the increased sensitivity of the cell voltage to the fraction of 
the hydrogen in the fuel channel. 


4.1.2. SOFC with fuel processed by CPOX 

Simulations are also performed for the case where the fuel 
supplied to the SOFC comes from a CPOX processor, and the 
operating parameters listed in the right column of Table 2 are 
considered. The approach used in [11] is adopted to calculate 
the composition of the reformate coming out of the CPOX as 
a function of temperature and oxygen-to-carbon ratio (yo,/c), 
where yo,/c is defined as follows: 


molar flow rate of O2 entering the reformer 


YO2/C = (37) 


molar flow rate of CH4 entering the reformer ` 

Fig. 4 illustrates simulation results of the steady-state cell 
voltage and solid structure temperature of the SOFC at differ- 
ent operating conditions. No significant difference between the 
baseline model and the MGFE model is observed, except for the 
condition with yfr and yo,/c as low as 0.2 and 0.5, respectively. 
Under these conditions, the cell voltage predicted by the MGFE 
model is about 0.02 and 0.04 V lower than the baseline model 
at current density of 0.3 and 0.6 Acm~’, respectively, and the 
solid structure temperature about 16 and 22 K lower, respec- 
tively. However, these operating conditions are not expected in 
real world operations, except in some short transient. For exam- 
ple, low Yo, /c can lead to carbon deposition in the reformer due 
to low H20/CH; ratio (see Fig. 5), and therefore yo,/c = 0.6 
is often suggested [12]. For high system efficiency, the fuel 
utilization is often kept above 0.8 [19]. 

Similar reasons as discussed for the first case attribute to 
the model error introduced by the MGFE approach in this 
case. Table 4 compares the gas composition in the fuel channel 
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Fig. 4. Steady-state voltage and temperature of SOFC in Case 2. 
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obtained by the baseline model and the MGFE model. One can 
see that Xcu, is also under-predicted by the MGFE model and 
the difference becomes more significant at lower yfr and yo,/c. 
The molar fractions of CH4 in the reformate fed into the fuel cell 
is shown in Fig. 5. As yo,/c decreases, more CH4 is unreformed 
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Fig. 5. Molar fraction of CH4 and H20/CH; ratio in the reformate of CPOX. 


in the CPOX and enters the fuel channel, leading to larger 
error in the composition calculation in the MGFE model. The 
over-counted SR reaction results in the reduced temperature in 
the SOFC and consequently the cell voltage, as shown in Fig. 4. 


4.2. Analysis of dynamic responses 


Dynamic responses of the two models are then compared 
and analyzed. Given that our primary interests lie in the con- 
trol system development and system integration for which the 
dynamic response of the system is very critical, we focus on 
the analysis of the model performance using both time-domain 
simulation and frequency domain analysis. The two cases with 
different fuel processing methods (i.e., SR and CPOX) are also 
considered here. 


4.2.1. Dynamic responses of SOFC models in Case 1 

Fig. 6 compares the open-loop responses of the two models 
in Case | to the step changes in the current load and gas supplies, 
based on simulation results. ig and we are initial operating cur- 
rent density and gas inlet flow rates, respectively, corresponding 
toi = 0.3 A cm~? and ye = 0.8. Two different pre-reforming 
ratios (Ypr = 0.2, 0.8) are considered. Up and Tọ are the steady- 
state operating voltage of the cell and the temperature in the solid 
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Table 4 
Gas composition in fuel channel of SOFC in Case 2 
i(Acm~*) YO2/C Ver Model Xcu, Xco, Xco Xpo Xm XN, 
0.8 Baseline 0.0005 0.1506 0.0393 0.3146 0.0652 0.4298 
0.6 MGFE 0 0.1518 0.0385 0.3124 0.0681 0.4293 
0.2 Baseline 0.0032 0.0510 0.1373 0.0977 0.2788 0.4321 
0.3 MGFE 0.0003 0.0517 0.1383 0.0932 0.2868 0.4296 
0.8 Baseline 0.0016 0.1612 0.0427 0.3337 0.0741 0.3866 
0.5 MGFE 0 0.1618 0.0431 0.3300 0.0798 0.3854 
0.2 Baseline 0.0107 0.0493 0.1493 0.0869 0.3102 0.3936 
MGFE 0.0014 0.0479 0.1561 0.0767 0.3314 0.3864 
0.8 Baseline 0.0009 0.1491 0.0406 0.3166 0.0627 0.4301 
0.6 MGFE 0 0.1502 0.0400 0.3139 0.0665 0.4293 
0.2 Baseline 0.0045 0.0496 0.1378 0.1008 0.2741 0.4332 
0.6 MGFE 0.0002 0.0500 0.1401 0.0948 0.2854 0.4295 
0.8 Baseline 0.0026 0.1600 0.0433 0.3368 0.0699 0.3874 
0.5 MGFE 0 0.1601 0.0448 0.3316 0.0781 0.3854 
0.2 Baseline 0.0145 0.0484 0.1480 0.0925 0.3002 0.3965 
MGFE 0.0008 0.0459 0.1585 0.0781 0.3307 0.3860 


structure at the initial operating setpoints, respectively. Despite 
the steady-state errors as discussed earlier, the dynamic trends 
exhibited in the two models are more consistent, as shown in 
Fig. 6, even for 7p; as low as 0.2. 

Frequency analysis are performed to compare the dynamics 
of the linearized baseline and MGFE models through the bode 
plots as shown in Fig. 7, where two operating conditions with 
different yp, are considered. Frequency responses of the solid 
structure temperature (To1) and the cell voltage (U) to the 
fuel inlet flow rate (Nin) and the air inlet flow rate (Nina) are 
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Fig. 6. Dynamic responses of the baseline and MGFE models in Case 1. 


considered. Here, Nin and Nin, are deemed as two primary 
inputs available for controlling the response of the SOFC. From 
Fig. 7, we conclude that: 


(1) Different yp, may lead to different values in gains and 
phases, but the correlation between the two models are 
similar at those two operating conditions. 

(2) The MGFE model has almost the same frequency response 
to Nin,a as the baseline model over a wide frequency range. 

(3) The frequency responses to Ning in the two models are 
matched well up to 0.1rads~!, while significant model 
errors appear at higher frequency. 

(4) The responses of 7,0) to the inputs may be approximated by 
first-order dynamics for frequencies lower than 0.1 rad s~!. 
The gains start to fall off at about 1073 rad s~!, which cor- 
responds to a time constant approximately equals to 1000 s. 


Since, as mentioned in Section 3, the dynamics of the air chan- 
nel flow in the MGFE model are kept the same as the baseline 
model, it is expected that these two models have quite matched 
frequency response to the air flow input Nina, as shown in 
Fig. 7. For the fuel flow, however, since the mass and energy 
balance dynamics in the gas flows are fast, the difference in the 
two models as mentioned earlier results in large model errors 
only at high frequency, as reflected in Fig. 7. This effect is more 
pronounced in the response of U, compared to Tsoi, because 
U can be directly affected by the gas composition in the fuel 
channel through the algebraic electrochemical relation. 

As an illustration, Fig. 8 compares the dynamic responses 
of the two models to a step increase in the fuel inlet flow rate, 
where both short and long scope responses are shown. The sim- 
ulation condition is the same as in Fig. 6 with ypr = 0.8. One 
can see that, while the responses of the molar fraction of H2 
and H20 (i.e., XH, and Xy,0) in the fuel channel exhibit fast 
dynamics in the baseline model, they change instantaneously in 
the MGFE model following the change in Nin s, resulting in the 
instantaneous change in the cell voltage. This difference is con- 
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Fig. 7. Bode plots of the linearized baseline and MGFE models in Case 1. 


sistent with the different dynamics between the two models at 
high frequency shown in Fig. 7. While the MGFE model exhibits 
matched behavior with the baseline model in long time horizon, 
it does not capture the very fast dynamics in the bulk flow in the 
fuel channel. 

This modeling error can be mitigated by introducing a ficti- 
tious filter, which has similar characteristics as the fast dynamics 
shown in Fig. 8, on the input of Nin. Fig. 9 compares the modi- 
fied MGFE model with the baseline model, where a filter with a 
time constant t = 0.5s is applied to the input of Nin f in the 
MGFE model. One can see that, with this modification, the 
MGFE model matches the baseline model much better at high 
frequency, compared to the results shown in Fig. 7. 


It is worth to point out that significant difference can also 
be found between high frequency responses of the MGFE and 
baseline models to the current (J) change which is considered 
as a disturbance input to the system in the future control design. 
Similar approach as mentioned above can be used to reduce 
the modeling error by modifying the MGFE model with simple 
fictitious filters. 

In practical SOFC systems, since the fuel flow has to pass the 
fuel supply and processing subsystems, including fuel reformer, 
heat exchangers and gas manifolds, before reaching the fuel 
channel of the SOFC stack, the dynamics associated with these 
components will act as the filter and therefore reduce or elim- 
inate the modeling error introduced by the MGFE assumption. 
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Besides, we are interested in the overall load following perfor- 
mance and the solid structure temperature in the SOFC, which 
are affected dominantly by the slow dynamics in Tso]. Therefore, 
given the purpose the low-order model is intended to serve, the 
MGFE model preserves the key dynamics in the baseline model 
and is amenable for model-based control design and analysis of 
SOFC systems. 

Further inspection of Fig. 7 could also lead to the conclusion 
that the two inlet flow rates (Ning and Nina) impart different 
dynamic control authorities on the system responses. While T3o) 
exhibits slow dynamics to inputs with similar time constants, 
the cell voltage responds differently. Two different time scales 
can be clearly identified: a slow one corresponding to the ther- 
mal inertia and a fast one due to the change in gas composition 
through mass balance dynamics. Note that Nin has larger gains 
on U compared to Nina, Suggesting it be the primary authority to 
control the voltage response of the SOFC. Both Ning and Nina 
can be used for thermal management in the cell. 


4.2.2. Dynamic responses of SOFC models in Case 2 

Fig. 10 compares open-loop responses of the SOFC baseline 
model and the MGFE model with fuel processed by a CPOX 
reformer (Case 2) when step changes are introduced in the cur- 
rent load and gas inlets. Initial operating conditions include 
i=0.3Acm7?, yr = 0.8 and YO,/C = 0.6/0.5. Similar to Case 
1, the dynamic responses of the simplified model match those of 
the baseline model. The bode plots of the linearized models in 
this case are also compared in Fig. 11. Similar conclusions, as 
elaborated in Case 1, can be drawn that the two models have close 
dynamic characteristics at frequency up to 0.1 rads~!. While 
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Fig. 10. Dynamic responses of the baseline and MGFE models in Case 2. 
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Table 5 
Comparison of computation efficiency 
Model Number of states in SOFC model CPU time (s) 
Case 1 Case 2 
Ypr = 0.8 Ypr = 0.2 Yoc = 0.6 Yoc = 0.5 
Baseline 11 3.14 3.26 3.80 3.24 
MGFE 6 2.65 2.56 251 2.48 


significant model errors can be found at higher frequency for the 
MGFE model, this error can be effectively reduced by introduc- 
ing a simple filter on Nin ¢ in the MGFE model. When integrated 
with dynamic models of fuel processing and gas supply sub- 
systems, we can eliminate this fictitious filter as the manifold 
filling and other dynamics associated with the gas supply path 
will have the same effects. 

In Table 5, computation efficiencies of the MGFE model and 
the baseline model are compared based on the simulation con- 
ditions for Figs. 6 and 10. The CPU time is the average of three 
simulations and all simulations are run on a laptop with Pen- 
tiumM 1.4GHz CPU and 512MB memory. From Table 5, it can 
be seen that the computation time is reduced 23.6% on aver- 
age using the MGFE model, compared to the baseline model. 
Meanwhile, the simplified model reduce the total number of 
states from 11 to 6, making it more suitable for model-based 
control design and system optimization. 


5. Conclusions 


In this paper, the MGFE method is adopted to simplify the cal- 
culation of the composition and consequently the mass balance 
dynamics of the bulk flow in the fuel channel of the SOFC. Using 
this approach, 5 out of 11 states can be removed from the base- 


line model. Both the steady-state and dynamic responses of the 
simplified model and the baseline model are compared through 
simulations and analysis. Different compositions of the inlet fuel 
flow and operating conditions are considered to investigate the 
impacts of using the MGFE method on the modeling perfor- 
mance. Analysis results show that the reduced-order model is 
able to preserve the key dynamics exhibited by the baseline 
model. Good match between the two models in steady-state 
response is also observed at realistic operating conditions. 
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